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Abstract 

If supersymmetry exists at the electroweak scale, then it should be discovered at 
the LHC. Determining masses, of supersymmetric particles however, is more difficult. 
In this paper, methods are discussed to determine combinations of masses and of 
branching ratios precisely from experimentally observable distributions. In many cases 
such measurements alone can greatly constrain the particular supersymmetric model 
and determine its parameters with an accuracy of a few percent. Most of the results 
shown correspond to one year of running at LHC at "low luminosity", 10 33 cm~ 2 s _1 . 
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1 Introduction 



If supersymmetry (SUSY) exists at the electroweak scale, then gluinos and squarks will be 
copiously produced in pairs at the LHC and will decay via cascades involving other SUSY 
particles to the lightest SUSY particle (LSP), the Xi- 111 most models the xi is stable, must 
be neutral and therefore escapes the detector. It should then be easy to observe deviations 
from the Standard Model such as an excess of events with multiple jets plus missing energy 
]f/T or with like-sign dileptons £ ± £ ± plus $t 0? H> 0- Determining SUSY masses is more 
difficult because each SUSY event contains two LSP's, and there are not enough kinematic 
constraints to determine the momenta of these. 

The strategy developed in Ref. M and in this paper involves three steps. First, we use a 
simple inclusive analysis to establish a deviation from the Standard Model. We select events 
with at least four jets and large missing energy and plot the distribution of 

M eS = p Tjl + p Tt2 + PT,3 + PT,4 + $T ■ 

Typically, this is dominated by Standard Model processes at low M c g- but is a factor of 
5-10 larger than the Standard Model prediction for large M e s. The value of M e s at which 
the signal exceeds that standard model backgrounds provides a first estimate of the SUSY 
masses. 

The second step in the strategy is to identify characteristic signatures of particles oc- 
curring near the end of the SUSY decay cascades and to use these as the starting point for 
further analysis. This is best explained by an example. Suppose that gluinos are slightly 
heavier than squarks and that x!i ~^ Xxh is kinematically allowed. Then one can have the 
following decay chain: 

9 + 9 

i i 

i i 
xl + q Xi + q 
I 

Xi+h 
i 

b + b 

Such an event typically contains two hard jets from the (Jl,r decays, two b jets from the h 
decay, large fir, and soft jets from the gluino decays and from gluon radiation. In this case, 
one can reconstruct h — > bb as a peak in the bb mass distribution and measure its mass H. 
Then the h can be combined with either of the two hard jets. The smaller of these two 
masses must be less than the squark mass and so has a sharp edge that measures a known 
function of the qi, x\i an d Xi masses. In many cases several such measurements can be 
made to determine several combinations of masses more or less precisely. 

Given actual data, the third step would be to make a global fit of a SUSY model to all 
available measurements, including both the precision measurements just described and more 
inclusive ones such as the jet, lepton, and 6-jet multiplicities and pr distributions. Such 
an analysis involves simulating large numbers of signal samples and is beyond the scope 
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Table 1: SUGRA parameters for the five LHC points. 



Point 


m 


mi/2 


A 


tan p 


sgn/z 




(GeV) 


(GeV) 


(GeV) 






1 


400 


400 





2.0 


+ 


2 


400 


400 





10.0 


+ 


3 


200 


100 





2.0 




4 


800 


200 





10.0 


+ 


5 


100 


300 


300 


2.1 


+ 



of this study. Instead, we try to determine the SUSY parameters using just the precision 
measurements of combinations of masses. In some cases this almost completely determines 
the SUSY model, while in others it significantly constrains it. 

What precision measurements can be made is very dependent on the SUSY model and so 
must be studied for specific SUSY parameters. The ATLAS and CMS Collaborations at the 
LHC have been considering five points in the minimal supergravity (SUGRA) model listed in 
Table [I] [|J. The SUGRA model || has the minimal SUSY particle content; universal scalar 
masses mo, gaugino masses mi/2, and trilinear breaking terms Aq at the grand unification 
scale; and radiative electroweak symmetry breaking driven by the large top quark mass. 
After electroweak breaking, the remaining parameters are a ratio of vacuum expectation 
values tan (3 at the weak scale and a sign, sgn fj, — ±1. We assume a default value for the top 
quark mass of 175 GeV and comment on the sensitivity to it below. While this model may 
not be the one that nature has chosen, we would like to emphasize that simulations can only 
be performed in the context of a consistent model. This is because many promising signals 
that might be clearly distinguished from Standard Model backgrounds in one channel, can be 
obscured by production and decays of other supersymmetric particles. We may not believe 
in this model, but the model that nature has chosen will be self-consistent. 

In the SUGRA model, if \2 ~^ xV 1 * s kinematically allowed, it has a substantial branching 
ratio and provides one starting point. This is the case for LHC Points 1, 2, and 5. If this 
decay is kinematically forbidden, then in many cases x\ ~ * Xi^ + ^~ can be observed. This is 
the case for LHC Points 3 and 4. The endpoint of the £ + £~ mass distribution provides a direct 
measure of M(x°) — M(xi), and opposite-sign, same-flavor dileptons can be used to identify 
events containing M(x\)- Other modes exploited in this paper include x\ ~^ ~~ * X?^ ± ^ =F 
and x\ ~^ XiW T — » Xi e± A t=Fz/ ^- 

All the analyses presented here are based on isajet 7.22 and a toy detector simulation. 
At least 50k events were generated for each signal point. The standard model background 
samples contained 250k events for each of ti, WZ with W — > eu, fiu, tu, and Zj with 
Z — > uu, tt, and 5000K QCD jets (including g, u, d, s, c, and b) divided among five bins 
covering 50 < pr < 2400 GeV. The calorimeter energy resolutions are taken to be 

EMCAL 10%/ VE + 1% 
HCAL 50%/y/E + 3% 
FCAL 100% /VE + 7%, H > 3 . 
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A uniform segmentation Arj = A0 = 0.1 is used with no transverse shower spreading. This 
is particularly unrealistic for the forward calorimeter. Jets are found using GETJET 
with a fixed cone size R = 0.4 or 0.7. Missing transverse energy is calculated by taking the 
magnitude of the vector sum of the transverse energy deposited in in the calorimeter cells. 
The jet multiplicity in SUSY events is rather large, so we will use a cone size of 0.4 unless 
otherwise stated. A lepton efficiency of 90% and a 6-tagging efficiency of 60% is assumed. 
Isolated leptons are required to satisfy an isolation requirement that no more than 10 GeV of 
additional Et be present in a cone of radius R = 0.2 around the lepton. Results are presented 
for an integrated luminosity of 10fb _1 , corresponding to one year of running at 10 33 cm _2 s _1 
so pileup has not been included. We will occasionally comment on the cases where the full 
luminosity of the LHC, i.e. 10 34 cm _2 s _1 , will be needed to complete the studies. For many 
of the histograms shown, a single event can give rise to more than one entry due to different 
possible combinations. When this occurs, all combinations are included. 

In Section |^ of this paper, we discuss using M eff to get a rough estimate of SUSY masses. 
We then turn to more detailed analyses. In these, we shall usually make cuts so that the 
Standard Model backgrounds are very small. These cuts are not optimal, particularly in 
the case of the higher mass points where event rates are lower. It may be desirable to have 
more signal events at the cost of more background. In Sections |3] and [| we discuss LHC 
Points 3 and 4 respectively. These points have rather light gauginos, so x\ ~^ Xi^ + ^~ can 
be used to measure M(x°) — M(xD- In Section || we discuss LHC Point 5, which has both 
X% ~ ¥ X\h ~~ * Xibb an d X2 ~^ — * Xi^ ^ ■ m Section | we briefly discuss LHC Points 1 
and 2. These have gluino and squark masses of about 1 TeV, so they really need more than 
10 fb" 1 for detailed study. In Section [7] we investigate how well the precision measurements 
discussed determine the parameters of SUSY model. Finally, we draw some conclusions. 
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2 Effective Mass Analysis 



The first step in the search for new physics is to discover a deviation from the Standard 
Model and to estimate the mass scale associated with it. SUSY production at the LHC is 
dominated by gluinos and squarks, which decay into multiple jets plus missing energy. A 
variable which is sensitive to inclusive gluino and squark decays is the effective mass M e s, 
defined as the scalar sum of the p^s of the four hardest jets and the missing transverse 
energy $ Tl 

M eS =p T ,l+ PT,2 + PT,3 + VT£ + $T ■ 

Here the jet pr's have been ordered such that pt,\ is the transverse momentum of the leading 
jet. The Standard Model backgrounds tend to have smaller $t, fewer jets and a lower jet 
multiplicity. In addition, since a major source of ft? is weak decays, large $t events in the 
Standard Model tend to have the missing energy balanced by leptons. To suppress these 
backgrounds, the following cuts were made: 

• $ T > 100 GeV; 

• > 4 jets with p T > 50 GeV and p TA > 100 GeV; 

• Transverse sphericity St > 0.2; 

• No fi or isolated e with px > 20 GeV and r\ < 2.5; 

• $ T > 0.2M eff . 

With these cuts and the idealized detector assumed here, the signal for all five LHC points 
is much larger than the Standard Model backgrounds for large M e g, as is illustrated in 
Figures. |l|-|5|. 

The peak of the M e s mass distribution, or alternatively the point at which the signal (S) 
begins to exceed the standard model background (B), provides a good first estimate of the 
SUSY mass scale, which is defined to be 

Msusy = min(M § , M fi J 

While M SU sY obviously should be set by the gluino and squark masses, the choice of Mu R 
as the typical squark mass is somewhat arbitrary. The ratio of the value M e fj for which 

Table 2: The value of M eff for which S = B compared to M S usy, the lighter of the gluino 
and squark masses. Note that Point 3 is strongly influenced by the $t and jet pt cuts. 



LHC Point 


M eff (GeV) 


Msusy (GeV) 


Ratio 


1 


1360 


926 


1.47 


2 


1420 


928 


1.53 


3 


470 


300 


1.58 


4 


980 


586 


1.67 


5 


980 


663 


1.48 
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S = B to Msusy was calculated by fitting smooth curves to the signal and background and 
is given in Table |[ It must be noted, however, that for LHC Point 3 the cuts produce a 
minimum value of M e g near the crossover. A more realistic treatment of the fa resolution 
could be important for this point. At this point event rates are so large that this step in our 
procedure is not needed; we will not use fa in the analyses shown below. 

To see whether the approximate constancy of this ratio might be an accident, 100 SUGRA 
models were chosen at random with 100 < m < 500 GeV, 100 < mm < 500 GeV, —500 < 
A < 500 GeV, 1.8 < tan/5 < 12, and sgn/i = ±1. These models were compared to 
the assumed signal, LHC Point 5. The light Higgs was assumed to be known, and all 
the comparison models were required to have Mh within ±3 GeV of its nominal value; the 
3 GeV error is determined not be the experimental measurements but by an estimate of the 
theoretical uncertainties on the prediction of Mh from the parameters of the SUGRA model. 
A sample of IK events was generated for each point, and the peak of the M ef j distribution 
was found by fitting a Gaussian near the peak. Figure |5] shows the resulting scatterplot of 
-^susy vs. M cS . The ratio is constant within about ±10%, as can be seen from Figure [?|. 
The mean value of the ratio is higher here than in Table |2| because this analysis uses the 
peak of the distribution rather than the point at which it is equal to the background. The 
error on the ratio is conservative, since there is considerable contribution to the scatter from 
the limited statistics and the rather crude manner in which the peak was found. While M e g 
does not provide a precise measurement, it has the advantage of being generally applicable 
to a broad range of SUSY models. 
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Figure 1: LHC Point 1 signal and Standard Model backgrounds. Open circles: SUSY signal. 
Solid circles: tt. Triangles: W — > £v, TV. Downward triangles: Z — > uu, tt. Squares: QCD 
jets. Histogram: sum of all backgrounds. 




Figure 2: SUSY signal and Standard Model backgrounds for LHC Point 2. See Figure 1 for 
symbols. 
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LHC Point 3 
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Figure 3: SUSY signal and Standard Model backgrounds for LHC Point 3. See Figure 1 for 
symbols. 
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Figure 4: SUSY signal and Standard Model backgrounds for LHC Point 4. See Figure 1 for 
symbols. 
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Figure 5: SUSY signal and Standard Model backgrounds for LHC Point 5. See Figure 1 for 
symbols. 
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Figure 6: Scatterplot of Msusy = min(Mg, M„) vs. M e s for randomly chosen SUGRA models 
having the same light Higgs mass within ±3 GeV as LHC Point 5. 
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Figure 7: Ratio M e ff/M SU sY from Figure []. The distribution of this ratio is approximately 
Gaussian with a width about 10% of its mean. 
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3 LHC Point 3: ra = 200 GeV, m 1/2 = 100 GeV, tan f3 = 2 



LHC Point 3 has relatively light superpartners and hence a very large production rate. At 
this point 200,000 events were generated, corresponding to about 1 week of LHC running at 
low luminosity; the statistical fluctuations on the plots are due to this small Monte Carlo 
sample. All of the squarks of the first two generations are heavier than the gluino, but 
one of the stop and sbottom mass eigenstates is lighter than the gluino while the other is 
heavier. The dominant SUSY production process is gg, and the dominant gluino decay mode 
is g — > bb. The lightest b is mainly b^ and so decays principally into bx% srnce X2 (Xi) is 
mainly W$ (B). Then X2 decays via virtual sleptons to Xi e+e ~ with a 16% branching ratio. 
SUSY events at this point are therefore dominated by final states involving 6-jets and pairs 
of opposite-sign, same-flavor leptons. Missing transverse energy is not used in the analysis 
at this point. 



3.1 Measurement of M 



X2 



Events are selected by requiring: 



• A pair of isolated leptons of opposite charge and the same flavor with pri > 10 GeV 
and \r)g\ < 2.5; 

• At least two jets tagged as b quarks and having p t > 15 GeV and \r)\ < 2; a tagging 
efficiency of 60% is assumed. 

The dilepton invariant mass distribution is shown in Figure |^. The dominant Standard 
Model background is ti production, which is quite small because it has smaller color factors 
and requires two leptonic decays. This background, as well as the combinatorial background 
from events with two x decays, can be eliminated by calculating the subtracted distribution: 



da 



dM 



da 



sub 



dM 



+ 



da 



dM 



da 



re" 



dM 



da 



dM 



This subtracted mass distribution has a sharp edge at Mg+g- = M^ 2 — M^, enabling this 
mass difference to be measured with great precision. In view of the enormous size of the 
event sample, the uncertainty on this measurement will be limited by systematic effects. The 
large sample of Z —>■ £ + £~ decays will be used for calibration both of the mass scale and 
of the relative e and \x acceptance. The methods employed will be similar to those used by 
CDF and DO in their determinations of the W mass 0, M. An estimate of 50 MeV for the 
uncertainty on Mj 2 — Mc 1 should be conservative. 



3.2 Gluino and Sbottom Reconstruction 

The next step is a reconstruction of the gluino and sbottom masses by combining a dilepton 
pair near the mass edge with jets. Events are selected that have 

• At least two jets, tagged as having a b quark with p t > 15 GeV and \rj\ < 2; a tagging 
efficiency of 60% is assumed; 
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• A e + e pair with 45 GeV < Mf+^- < 55 GeV and no other electrons or a pair 
in the same mass range and no other muons in the event. 

Since the mass of the lepton pair is near its maximum value, in the rest frame of X2 both xi 
and the £ + £~ pair are forced to be at rest. The momentum of X2 m the laboratory frame is 
then determined to be 

P i2 = (l + M^/M t+e -)P e+e -- 

where M^o must be assumed (see below). This momentum can be combined with a 6-jet to 
determine m~ b and a second 6-jet to determine m g . The 6-jet energy and momentum must 
be corrected for the fact that particles are lost outside the R = 0.4 jet cone and for the 
fact that weak decays produce neutrinos in the jets. In this study, the correction factor was 
determined using the data generated for LHC Point 5, where the Higgs peak (h — > 66) is 
observable. In practice, techniques similar to those of references [ID, |TT| would be used at 
LHC. 

Figure ^ shows a scatterplot of m g — m~ b vs. m g . Projections onto the axes, shown in 
Figures |10j and |TX , have clear peaks. The positions of the peaks determine m g — m~ b and m~ b 



assuming that M^o is known. Again, statistical errors are small and the dominant errors will 
be from the determination of the jet energy scale. A careful jet energy calibration has not 



been performed, so the peaks in Figures. [10] and [IT] are displaced slightly from their nominal 



values of 277.8 and 20.3 GeV. These systematic errors can be estimated from those currently 



obtained by CDF and DO in the determination of the top quark mass [[10j, |TT|. The mass 
difference m g — m b is insensitive to the assumed Xi mass while the reconstructed sbottom 
peak moves. 



The dependence of the b mass peak on the assumed value of M^o is shown in Figure |12| , 
where M^o is varied by ±20 GeV from its nominal value. In making this plot we have 
required that the mass difference M(x2&6) ~ ^(x^b) be within 15 GeV of the value where its 
distribution peaks. This cut removes considerable background as can be seen by comparing 
the peaks in this figure with that in Fig|lO|. We estimate 

M^(measured) — Mg(true) = 1.5 [M^o (assumed) — M^o(true)) ± 3 GeV 

and 

Mg(measured) - M^(measured) = M^(true) - Mg(true) ± 2 GeV 

The Xi mass will be determined by a global fit of the SUSY model to all the measurements; 
see Section 0. 



3.3 Light Squark Reconstruction 

Light squarks can also be reconstructed at this point using the decay chain qi — > x°g, which 
has a branching ratio of approximately 10%. There is an enormous background from gluino 
decays to bb, so events must be rejected if there is a 6-jet present. We use the ATLAS 
6-tagging study (see Figure 3.42 of Ref. |]]). At low luminosity this study implies that a 
tagging efficiency of 90% for 6-jets can be achieved at the price of misidentifying 25% of the 
light quark jets as 6-jets. While this mistag rate is not adequate in the cases where a 6-tag 
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is required, it implies that 90% of the 5-quark jets can be vetoed and 75% of the light quark 
jets accepted by the same cut. This veto prescription is used in this subsection. 
Events are selected as follows: 

• At least one jet with p t > 125 GeV and \rj\ < 2. 

• No 6-jets with p t > 15 GeV and \rj\ < 2; a vetoing efficiency of 90% is assumed and 
25% of non 6-jets are assumed to be rejected also. 

• An e + e~ pair with 45 GeV < M^i- < 55 GeV and no other electrons or a yU + /i~ pair 
in the same mass range and no other muons in the event. 

The reconstruction of the momentum of x® * s performed using the same method as above 
by selecting events near the endpoint of the dilepton mass distribution. We assume that 
the SUGRA model is used to infer the mass of Xi from the X2 ~ Xi mass difference. Jets 
of \r]\ < 2 and pt > 125 GeV are now combined with the X2 an d the mass distribution is 
shown in Figure |13[ Even with the 90% vetoing efficiency for b— quarks there are a significant 
number of b— jets remaining in this plot. The contribution from the light squarks is shown 
as the dashed-histogram. If the vetoing efficiency were raised to 95% approximately one-half 
of the remaining b— jets are removed and consequently the peak moves to a larger mass. The 
peak shown has contributions from of mass 278 GeV and the light quarks that have mass 
around 310 GeV. Charge —1/3 and +2/3 squarks are separated by about 5 GeV in mass; this 
contributes to the broadening of the peak. That the peak is real can be seen be estimating 
the combinatorial background as follows. Events are mixed by taking the X2 momentum 
from one event and the jet from another; both events satisfying the same selection criteria. 



The mass distribution obtained in this way is shown as the hatched distribution in Figure 13. 



Conservatively, we estimate an error of 20 GeV on the average qi mass from this method. 



3.4 Branching ratio of %2 — ► X\^ + ^ 

By selecting events with four tagged b— jets and either two or four isolated leptons, the 
product of branching ratios BR{x% ~^ Xi^ + ^~) x BR(b — > bx^X) can be determined. There 
are 150000 events/10 fb _1 with two dilepton pairs and four b— jets. The backgrounds from 
non supersymmetric sources are negligible, and again therefore the dominant uncertainties 
are systematic. Using a value of 3% for the uncertainty on the absolute lepton acceptance, 
we expect that BR(x2 X\ e+e ~) x BR(b — > bx%X) can be determined to be (14.0 ± 0.5)% 



3.5 Electroweak Production of Superpartners 

At this SUGRA point, sleptons cannot be produced from the decay of strongly interacting 
sparticles. The production rates are therefore quite small despite the low masses (m^ = 215 
GeV, rrie R = 206 GeV) as they must be pair produced in Drell-Yan like processes. The 
heavier charginos and neutralinos are only rarely produced in the decays of gluinos, so again 
their dominant production mechanism is electroweak. Unlike the case of sleptons, the direct 
production rate of the lighter charginos and neutralinos is quite large. An attempt has 
been made to isolate these processes. This is an example of a case where the analysis of a 
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complete SUSY signal is needed. The signals that we are attempting to extract stand clearly 
above Standard Model backgrounds, but we face the large background from the production 
of strongly interacting sparticles. As so few events pass the cuts, we generated separate data 
samples corresponding to the electroweak production of sparticles and reweighted the events 
appropriately. 

Events are selected that have: 

• A three isolated leptons a pair of which have opposite charge and the same flavor with 
p Tl > 10 GeV and < 2.5; 

• No jets with p t > 30 GeV in \r]\ < 3.0. 

The jet veto is needed to remove gluino and squark initiated events. These events have jets in 
the central region arising from the decay products of the sparticles and from final state gluon 
radiation. These events also have jets, approximately uniform in rapidity, from initial state 
radiation. This latter source is also present in the direct production of charginos, neutralinos 



and sleptons. Figure [14] shows the dilepton invariant mass distribution of the two leptons 
that have opposite charge and the same flavor. The number of generated events in this plot 
is not large, but are sufficient to demonstrate that in 10 fb _1 of data there will be sufficient 
events for a precise measurement. The background events in this plot (corresponding to 
three generated events) are from it production, the third lepton being from the decay of a 
6-quark. A stricter jet veto (20 GeV instead of 30 GeV) reduces this background further. 

There is an indication of an edge in the mass distribution corresponding to the decay 
X2 ~^ Xi^ + ^~ • The events in this plot are dominated by the production of x^Xi tinal states 
whose contribution is shown as the dotted histogram. If two isolated leptons are required 
and the same plot made the result is that there are more events. There is now a potential 
background from Drell-Yan production of dilepton events which must be eliminated by a cut 
on missing Ifix or the angle between the two leptons; the Drell-Yan events are back-to-back 
while in the SUSY events the leptons arise from \X ~^ Xi^ + ^~ an d are therefore close in 
angle. The production rates in these two and three lepton final states can be compared 
and used to provide a powerful argument concerning the origin of the lepton samples and 
provide an additional constraint on the model since, as we will demonstrate in section [7], the 
measurements that have been made using the strong production of sparticles fix the model 
parameters, resulting in a prediction for the rates shown in Figure [TJ]. 

In principle, the decay ei — > x\ e should be reconstructible by selecting with a least 3 
isolated leptons, an oppositely charged pair of which have mass between 45 and 55 GeV. 
The momentum of x° is reconstructed as above and then combined with a third lepton to 
search for a reconstructed e^. The extraction of this signal is very difficult. The production 
rate for gauginos provides a serious background that can only be controlled by increasing 
the number of isolated leptons required. The dominant slepton production process is e~L + v e - 
This can be extracted only by requiring at least jour isolated leptons from the decay chain 

it + h 

I I 

i i 
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or alternatively from 

n + vi 
i i 

I I 

X\ + i + + I- Xi + ^ + ^ 

The dominant decay chain z>£ — > x^f, Xi" ~~ *■ X? + jets is killed by the jet veto requirement. 
The experiment is only feasible at high luminosity. 
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Figure 8: The invariant mass distribution of e + e and jj, + fi pairs arising at Point 3. The 
background, shown as a hatched histogram is mainly due to tt events. 
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Figure 9: The reconstruction of gluino and sbottom decays from the decay chain g —> X2{~^ 
Xi£ + £~)b. Events are selected near the endpoint of the £~£ + mass distribution (mass between 
45 and 55 GeV) and the momentum of \2 reconstructed. Two fo-jets are then required and 
the mass of b + X2 (in = m~ b ) and the mass difference 5m = m bb ^ 2 — m b ^ 2 is computed. The 
scatterplot in these two variables is shown. The tVjet energies have been recalibrated and a 
tagging efficiency of 60% per b included. 
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Figure 10: The M(b) projection of Figure 
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Figure 11: The M(g) — M(b) projection of Figure |9|. The dashed histogram shows the 
projection if a cut is made requiring that the events lie in a slice of on the abscissa of 
between 230 and 330 GeV of Figure 0. 
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LHC Point 3 




Figure 12: The same as Figure 10 with the addition of two more histograms (dashed and 
dotted) showing the result if the assumed value of m^o is varied by ±20 GeV. A cut is 
imposed on the mass difference \rrig — m~ b — 20 1 GeV < 15 GeV before the projection of the 
scatterplots is made. 




Figure 13: Reconstructed q~i mass at Point 3. The combinatorial background estimate is 
shown as a hatched histogram and the events due to light squarks as the dashed histogram. 
The remaining events are due to gluino decays where a 6-jet is misidentified as a light quark 
jet. 
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Figure 14: The invariant mass distribution of e + e~ and pairs arising at Point 4. Events 
are selected requiring no jets with p t > 30 GeV in \r)\ < 3 and at least three isolated leptons, 
two of which are of the same flavor and opposite charge. Lepton detection efficiency of 90% 
per lepton is included. The dashed histogram shows the contribution arising from the direct 
production of X1X2 fi na l states. The background is shown as the hatched histogram. 
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4 LHC Point 4: ra = 800 GeV, m 1/2 = 200 GeV, tan f3 = 10 



LHC Point 4 has squarks which are much heavier than gluinos, so production of the latter is 
dominant. The heavier chargino and neutralinos have a much larger admixture of gauginos 
than at the other points. Hence, the gluinos decay into all combinations of all the charginos 
and neutralinos and all the quark pairs with comparable branching ratios, giving a very 
complex mixture of signatures. 



4.1 Selection of gaugino decays Xi ~^ Xj^ + £ 

The objective of this analysis is to isolate opposite-sign, same-flavor dileptons coming from 
X2 ~^ X34 — > XiZ, and xt ~^ xfZ- First, the following cuts were made to suppress 

the Standard Model backgrounds: 

• M eS > 800 GeV; 

• $ T > max(100GeV,0.15M eff ); 

• > 4 jets with p TA > 100 GeV, p T ,2,3,4 > 50 GeV; 

• £ + £~ pair with p T<e > 10 GeV, rj £ < 2.5; 

• i isolation cut: E T < 10 GeV in R = 0.2; 

• Transverse sphericity St > 0.2. 

The opposite-sign, same-flavor and opposite-flavor dilepton mass spectra for the Point 4 
signal and the opposite-sign, same-flavor Standard Model background with these cuts are 



shown in Figure [15|. There are clear low-mass and Z opposite-sign, same- flavor signals. (Note 
that the Z is treated as a narrow resonance in the event generator.) The same-sign Standard 
Model background is not shown but is smaller than the opposite-sign background. 

The difference of the opposite-sign, same-flavor and opposite-sign, opposite-flavor dilep- 
ton distributions is shown in Figure |16[ This difference should only have contributions from 
X2 ~^ Xi^ + ^~ (shown as a dashed curve) and from Z — £ + £~ decays from heavy charginos and 
neutralinos, which contribute to the Z peak. Contributions from two independent chargino, 
top, or W decays contribute equally to both flavor combinations and therefore cancel in this 
figure. The Standard Model background in the figure fluctuates in sign because of limited 
statistics but should also mostly cancel. 

The edge of the x° > Xi^ + ^~ signal is not quite as sharp as in previous cases, but it 
clearly can be measured with an error of ~ 1 GeV or less. The observation both of this edge 
and of the Z peak shows that both light and heavy gauginos contribute, since for any given 
gaugino Xh decay into XjZ is much larger than decay into Xj^ + ^~ ■ 

The relative number of events in the two parts of the distribution can be measured with a 
statistical error of a few percent. The systematic error on the e and [l acceptance should be 
comparable to or less than this after Z decays are studied carefully and used for calibration 
of the calorimeter. Since the sleptons are also heavy at this point, the leptonic branching 
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ratios for the Xi are essentially determined by the Z branching ratios, so the relative number 
of events provides a measure of 

E _ ± B(g - XiXMiCi - ZX) 



B(g - X2X) 

There are of course non-negligible corrections from squark production and from lepton ac- 
ceptance. 

To see how useful such a branching ratio measurement might be, samples of 10K events 
each with m = 800, 700, and 600 GeV were generated, forcing the decays X2 Xi& + £~ 
and Z — > . Figure [L7| shows the resulting opposite-sign, same- flavor mass distributions 
including the 6% leptonic branching ratios. (Subtraction of the same-sign background does 
not work properly when decays are forced, so the subtracted distribution is not shown.) The 
number of events in the Z peak and below 70 GeV are 210/2000, 260/2350, and 3000/3150 
respectively. The ratio is nearly constant, and the change in absolute number is similar to 
the change in the total cross section. There is, however, sensitivity to tan/?, which is more 
directly related to the mixing of the heavy and light gauginos. A sample of 10K events with 
mo = 800 GeV and tan (3 — 5 was generated, again forcing the decays. Figure shows the 
two distributions; the ratio is 80/1950. More study is needed, but it seems likely that this 
ratio could constrain tan (3 — 10 to ~ 10%. 

The results at this point are very sensitive to the top quark mass. For example, Figure |i~9] 



shows the same distribution as Figure 16 for the same SUSY parameters but for m t = 
170 GeV instead of 175 GeV. Note that the Z peak is dramatically larger. The reason 
for this extreme sensitivity is that this point is very close to the boundary of the allowed 
region: there is no electroweak symmetry breaking for m t = 165 GeV. The possibility of such 
sensitivity, however, suggests that theoretical uncertainties could play an important role in 
this region of parameter space. 

4.2 Selection of xl xfW T -> e ± /i T X 

Isolated e/x events come from two independent W or wino decays, not from single Z or 
neutralino decays. At LHC Point 4 the branching ratio for X_l -> XiW T is about 84%; this 
decay contributes to e ± i^ but not to e ± \i ± . Since the gluino is a Majorana fermion, other 
channels involving two independent gg or gq decays contribute equally to e ± n^ and e ± n ± . 

There is a large background to e ± fjL :f from ti production. To suppress this it is necessary 
to raise the cut on M eff from 800 GeV to 1000 GeV. The like-sign and opposite-sign e/x mass 
distributions with this cut and the other cuts described in the previous subsection are shown 
in Figure EH together with the Standard Model backgrounds. The difference of the like-sign 



and opposite-sign distributions is shown in Figure ^1| The Standard Model background 
shows statistical fluctuations but is fairly small after these cuts. 

For Point 4 the endpoint of the e ± fi T mass distribution is determined by x°a — > xfW T — ► 
e ± fi T Xi an d is 220.6 GeV. This is consistent with Figure ^T|. Of course other models might 
lead to the dominance of this decay by other modes. Because there are two missing neutrinos 
in addition to the Xii there is no sharp edge at the kinematic limit, so the endpoint can be 
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determined only roughly. However, the total number of events in Figure ^ can be measured 
to a few percent; it provides a measure of another combination of branching ratios. 
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Figure 15: M(+i- distribution for opposite-sign, same-flavor dileptons for the Point 4 signal 
(solid histogram), opposite-sign, opposite- flavor dileptons (dashed histogram), and Standard 
Model background (shaded histogram). 
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Figure 16: Difference of the M e + e - distribution for opposite-sign, same- flavor dileptons and 
opposite-sign, opposite-flavor dileptons for the Point 4 signal (open histogram) and the 
Standard Model background (shaded histogram). 
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Figure 17: Comparison of M^+£- distributions for opposite-sign, same-flavor dileptons for 
m = 800 GeV (solid histogram), m = 700 GeV (dashed histogram), and m = 600 GeV 
(dotted histogram). 
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Figure 18: Comparison of M^+£- distributions for opposite-sign, same- flavor dileptons for 
tan/3 = 10 (solid histogram) and tan/3 = 5 (dashed histogram). 
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LHC Point 4, m t = 170 GeV 
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Figure 19: The same as Figure [TBI but with m t = 170 GeV. Note the large change in the size 
of the Z peak. 
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Figure 20: M e ±^ (solid) and M e ±^± (dotted) distributions for the Point 4 signal, and Stan- 
dard Model opposite-sign, opposite- flavor background (shaded histogram). 
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Figure 21: Difference of the M^+£- distribution for opposite-sign, opposite-flavor and same- 
sign, opposite-flavor dileptons for the Point 4 signal (open histogram) and the Standard 
Model background (shaded histogram). 
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5 LHC Point 5: ra = 100 GeV, m 1/2 = 300 GeV, tan/3 = 2.1 



LHC Point 5 has a gluino with mass 767 GeV and light squarks with masses of 662- 
690 GeV, so g -> qq dominates. It has M(xjj) = 232.6 GeV, M(x°i) = 121.7 GeV, and 
M(h) = 104.15 GeV, so x\ xV 1 is kinematically allowed. It also has light right-handed 
sleptons, M(£r) = 157.2 GeV, so that x° ~~ > ^r^ t ~^ Xi^ + ^~ a ^ so h as a large branching ratio. 
This point was chosen so that the Xi provides the correct amount of cold dark matter for 



cosmology; this generally requires relatively light sleptons [12 



5.1 Selection of h — > bb and Measurement of M(u£) — M(xT) 

For LHC Point 5 the decay chain x\ Xi^, h — > bb has a large branching ratio, as is 
typical if this decay is kinematically allowed. The decay h — ► bb thus provides a handle for 
identifying events containing x^'s 0- Furthermore, the gluino is heavier than the squarks 
and so decays into them. The strategy for this analysis is to select events in which one 
squark decays via 

Q ~> xll, X°2 ~> Xih, h^bb, 

and the other via 

q^xlq, 

giving two b jets and exactly two additional hard jets. 

IS A JET 7.22 [0] was used to generate a sample of 100K events for Point 5, corresponding 
to about 5.6 ftT 1 so the signal statistics shown in this section roughly correspond to the 
actual statistics expected in 1 year at low luminosity. The background samples generally 
represent a small fraction of an LHC year. The detector response was simulated using the 
toy calorimeter described above. Jets were found using a fixed cone algorithm with R = 0.4. 
The following cuts were imposed: 

• $ T > 100 GeV; 

• > 4 jets with p T > 50 GeV and p T>1 > 100 GeV; 

• Transverse sphericity St > 0.2; 

• M eS > 800 GeV; 

• $ T > 0.2M eff . 

As before, jets were tagged as 6's if they contained a B hadron with p T > 5 GeV and r\ < 2, 
and a tagging efficiency of 60% per b was included. 

Figure |22] shows the resulting bb mass distributions for the signal and the sum of all 
Standard Model backgrounds with pT,b > 25 GeV together with a Gaussian plus quadratic 
fit to the signal. The mistagging background is comparable to the real background shown. 
The energy calorimeter scale for 6-jets was recalibrated to bring the Higgs mass peak to it's 
correct value which will be measured ultimately via the decay to 77. The correction is about 
8%. Using a larger cone, R = 0.7, gives an uncorrected peak which is closer to the true mass 
but wider. Note that for Point 5 the light Higgs could be discovered in this channel with 
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much less integrated luminosity than is needed to observe h — *> 77; the latter would provide 
a better mass measurement, AMh < 1 GeV. 

Events were then required to have exactly one bb pair with invariant mass within ±1.5(7 
(~ 19 GeV) of the Higgs peak and exactly two additional jets with p T > 75 GeV. The 
invariant mass of each jet with the bb pair was calculated. For the desired decay chain, one 
of these two must come from the decay of a single squark, so the smaller of them must be 
less than the kinematic limit for the decay chain q — > X2I ~* Xi^q, namely 



max N 2 
hq ) 



Ml + (M ( 



Mi* + Mj 



M k + 



(M| 



Ml 



2 \2 



AMIM^ 



2M? 



Using the average of the ul and di masses gives M^ ax = 506 GeV. The smaller of the two 
bbj masses is plotted in Figure ^ for the signal and for the sum of all backgrounds and 
shows an edge near the expected value. The Standard Model background shows fluctuations 
from the limited Monte Carlo statistics but seems to be small near the edge, at least for the 
idealized detector considered here. There is some background from the SUSY events above 
the edge, presumably from other decay modes and/or initial state radiation. 

A detailed understanding of the shape of this edge and its relation to the masses involved 



requires more study. Based on the statistics in Figure |23|, it seems likely that one could 
determine the end point of the spectrum to ~ 40 GeV in one year and to half of that in three 
years at low luminosity. 



5.2 Selection of W — > qq and Measurement of M{ul) — M(xi] 
Point 5 also has a large combined branching ratio for one gluino to decay via 

9 -> q L q, q L -> xf g, xf -> X?W /± , W ± 

and the other via 



qq- 



9 -> qaq, qn 



xlq, 



x°iq 



giving two hard jets and two softer jets from the W. The branching ratio for qi 
is small for Point 5, so the contributions from g —>■ qLq and from qL<iL pair production are 
suppressed. 

The same signal sample was used as in the previous subsection. The combinatorial 
background for this decay chain is much larger than for the previous one, so harder cuts are 
needed: 

• $ T > 100 GeV; 

• > 4 jets with p Tlj2 > 200 GeV, p T 3,4 > 50 GeV, and r] 3A < 2; 

• Transverse sphericity St > 0.2; 

• M eS > 800 GeV; 

• $ T > 0.2M eff . 
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The same 6-tagging algorithm was applied to tag the third and fourth jets as not being b 
jets. In practice one would measure the 6-jet distributions and subtract them. 

The mass distribution M 34 of the third and fourth highest px jets with these cuts is shown 



in Figure ^4] for the signal and the sum of all backgrounds. A peak is seen a bit below the W 
mass with a fitted width smaller than that for the h in Figure ^; note that the W natural 
width has been neglected in the simulation of the decaysjj] The Standard Model background 
is more significant here than for the h — > bb channel. Events from this peak can be combined 
with another jet as was done for h — > bb, providing another determination of the squark 



mass, Figure |25], with an error similar to the previous one. Figure also provides a starting 
point for measuring the W production rate in SUSY events. Knowing this rate is essential 
when searching for excess leptons from other sources such as gaugino decays. 

5.3 Selection of x° 2 ^ U ^ x\ii 

Point 5 has relatively light sleptons, as is generically necessary if the x? is to provide ac- 
ceptable cold dark matter fl2| , since Xi ~ Xi annihilation in the early universe proceeds via 
slepton exchange and the slepton mass must therefore be small enough to make this rate 
sufficiently large. Hence the two-body decay 

xl - - xTr 

is kinematically allowed and competes with the x® X\h decay, producing opposite-sign, 
like-flavor dileptons. This source of sleptons is much larger than that from direct production 
and their discovery is much easier than at Point 3 despite the fact that the total sparticle 
production rate is much larger at that point. The largest Standard Model background is tt. 
To suppress this and other Standard Model backgrounds the following cuts were made on 
the same signal and Standard Model background samples used previously: 

• M eS > 800 GeV; 

• $ T > 0.2M eff ; 

• > 1 R = 0.4 jet with p T]1 > 100 GeV; 

• £ + £~ pair with p T<e > 10 GeV, m < 2.5; 

• I isolation cut: E T < 10 GeV in R = 0.2; 

• Transverse sphericity St > 0.2. 

With these cuts very little Standard Model background survives, and the M u mass distri- 



bution shown in Figure ^6] has an edge near the kinematic limit for this decay sequence, 
namely 



1 £_ 1 2 

\ M\ \ M 



1 - — w 108.6 GeV 
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Table 3: Comparison of masses relevant to dilepton spectrum for Point 5 and for modified 
point with m = 120 GeV. 



Mass 


Point 5 


With m = 120 GeV 


9 


767.1 GeV 


767.2 GeV 




231.2 GeV 


231.4 GeV 


Ir 


157.2 GeV 


170.6 GeV 


xl 


121.3 GeV 


121.4 GeV 



Observing both h —>■ bb with Mh > Mz and an £ + £~ continuum with Mu > Mz would 
certainly suggest, and perhaps establish, the existence of light sleptons. 

If Mu is near its kinematic limit, then the velocity difference between the £ + £~ pair and 
the Xi is minimized in the rest frame of Having both leptons hard requires M^/M? ~ 
M^o/Mji. Assuming this and M^o = 2M^o implies that the endpoint in Figure ^ is equal 
to the Xi mass. An improved estimate could be made by detailed fitting of all the kine- 
matic distributions. Events were selected with M^ x - 10 GeV < M u < M^\ and the xl 
momentum was calculated using this crude Xi mass and 

p^a = (M^oJ M u ) p u ■ 

The invariant mass M u ^ of the £ + £~, the highest pt jet, and the Xi was then calculated 
and is shown in Figure |27|. A broad peak is seen near the light squark masses, 660-688 GeV. 
This peak contains complementary information to that obtained from Figure ^ and more 
information about the masses could be obtained by performing a combined fit. 

A detailed analysis of the dilepton mass spectrum would require varying M{xX)i M(£r), 
M(xi), and the Pt{xX) distribution and fitting the distributions of M(££), pr(££), Pr(h), and 
Pt{£'i) Ivt{£i)- Rather than do this, a sample of 50K events with mo = 120 GeV but otherwise 
the same parameters was generated. A comparison of some of the relevant masses is shown 
in Table the slepton mass changes by 13 GeV, and the rest are essentially identical. The 
M(££) mass distribution near the edge is shown in Figure There is a shift in the location 
of the edge by about 2 GeV, which should be observable. There is not much change in the 
Pt distribution of the ££ pair, Figure [2D|. The most sensitive distribution is Pt{^2)/pt(^i), 
where by definition pt{£?) < Pt{£i)- There is a clear change in the shape, as one would 
expect. For fixed values of the other parameters, one ought to be able to determine mo to 
~ 5 GeV, although careful subtraction of Standard Model backgrounds will be necessary. 

5.4 Top Production in SUSY Events 

Gluino decays at LHC Point 5 have a sizable branching ratio to ti. We have attempted to 
isolate such a sample by searching for top decays. In order to reduce the background from 
Standard Model top production, the cut on M e g has been raised to M e g > 1000 GeV for this 
analysis. The event selection is as follows: 

^The Higgs width is much smaller than the W width for Higgs masses relevant to these analyses. 
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• $ T > 100 GeV; 



• > 4 jets with p TX > 100 GeV, Pt2,3,4 > 50 GeV; 

• Transverse sphericity St > 0.2; 

• M cS > 1000 GeV; 

• $ T > 0.2M eff . 

• Exactly two b-jets with pt > 25 GeV. 

If the invariant mass of the two- 6 system is within ±20 GeV of the Higgs mass, the event is 
rejected. To search for W —>■ qq, we calculate the invariant mass of all dijet combinations with 
both jets having pt > 50 GeV and with neither jet being tagged as a b. Combinations with 
\M q q — M w \ < 10 GeV mass are considered W candidates. Combinations where 20 GeV < 
\M q q — M w \ < 30 GeV are used to model the shape of the background; we will refer to these 



events as the W sideband region. The M qq - distribution is shown in Figure |31] where a clear 
W peak is visible. The average jet multiplicity in these events is ~ 9 and hence there is a 
large combinatorial background in this figure. 

Dijets in the W mass and W sideband regions are then combined with each of the two 
b jets, giving two combinations for each W. Results of this reconstruction for the W signal 
and sideband regions are shown in Figure [I2| Subtraction of the sideband distribution from 
the signal yields Figure [33|. A clear top signal is seen. There is a small background of 
Standard Model production of top quarks, shown as the hatched histogram on this figure. 
The dominant decay chain giving rise to this signal is 

9 -> tt 
t - xfb 
xt - X°iW 

The kinematics of the decay t — > xfb, Xi ~* i&V restricts the Wb invariant mass distribu- 
tion to be between 203 and 356 GeV. There is some evidence in Figure ^3| for an excess of 
events in this region. Unfortunately, other possible combinations and decays such as g — > bb, 
b — > xtt> xt ~ * Xi^ 7 ) which has a smaller combined branching fraction, make the exact 
interpretation difficult. Thus, when discussing the determination of SUSY parameters in 



Section 7.2, we assume that top decays in SUSY events can be seen at Point 5 (indicating 



that the g —> tt channel is open), but we will not assume that the t mass can be determined. 
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Figure 22: M(bb) for pairs of b jets for the Point 5 signal (open histogram) and for the sum 
of all backgrounds (shaded histogram) after cuts described in the text. The smooth curve is 
a Gaussian plus quadratic fit to the signal. The light Higgs mass is 104.15 GeV. 
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Figure 23: The smaller of the two bbj masses for the signal and background events with 
73 < M(bb) < 111 GeV in Figure ^ and exactly two additional jets j with p T > 75 GeV. 
The endpoint of this distribution should be approximately M^ ax = 506 GeV. 
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Figure 24: M 34 for non-6 jets in events with two 200 GeV jets and two 50 GeV jets for the 
Point 5 signal (open histogram) and the sum of all backgrounds (shaded histogram). 
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Figure 25: The smaller of the two qqj masses for signal and background (shaded) events with 
71 < M(bb) < 87 GeV in the previous figure and with exactly two additional jets j with 
Pt > 75 GeV. The endpoint of this distribution should be approximately the mass difference 
between the squark and the Xi, about 565 GeV. 
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Figure 26: M« for the Point 5 signal (open histogram) and the sum of all backgrounds 
(shaded histogram). 



80 



v 60 
o 

> 

CD 

CD 

o 40 

£J 
^5 



LHC Point 5 



~i — i — i — | — i — i — i — | — i — i — r 



w 

P1 

P2 
P3 
P4 
P5 



df 11.50/9 
32.18 
623.5 
64.11 
-24.03 
.1042. 
-.7275E-04 



A 



20 



I I I I I ll I I I I I I I I I I I I I L_ 

200 400 600 800 1000 




M lljx (GeV) 



Figure 27: M U j^o for events with 86 < M u < 109 GeV using p^o = M^/M u f> u for the 
Point 5 signal (open histogram) and the Standard Model background (shaded histogram). 
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LHC Point 5, m = 1 00, 1 20 GeV 
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Figure 28: M(££) for Point 5 (solid curve) and for a modified point with m = 120 GeV 
(dashed curve). The shaded histogram is the sum of all Standard Model backgrounds. 




Figure 29: pr(^) for Point 5 (solid curve) and for a modified point with m = 120 GeV 
(dashed curve). The shaded histogram is the sum of all Standard Model backgrounds. 
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LHC Point 5, m = 100, 120 GeV 
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Figure 30: Pt{^2) /pr(^i) for Point 5 (solid curve) and for a modified point with m = 120 GeV 
(dashed curve). The shaded histogram shows the sum of all backgrounds. 




Figure 31: The invariant mass distribution used to search for W —> qq decays in the Point 5 
top reconstruction. The shaded histogram shows the sum of all backgrounds. 
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LHC Point 5 




Figure 32: The invariant mass distribution of qqb candidates for combinations where the qq 
pair are in the W mass region (solid) and W sideband region (dashed). The 6-jet energies 
have not been recalibrated and a tagging efficiency of 60% per b included. Note that each 
event appears twice. 
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Figure 33: The sideband subtracted invariant mass distribution of Wb candidates. The 
hatched distribution shows the contribution from Standard Model production of top quarks. 
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6 LHC Points 1, 2: ra = m 1/2 = 400 GeV, tan (3 = 2, 10 

LHC Points 1 and 2 have gluino masses of about 1 TeV and squark masses about 50 GeV 
lighter. These are close to the upper limit of the expected range if SUSY is to be relevant 
to electroweak symmetry breaking. The cross sections are quite small, so more than 10 ftT 1 
is needed for precision studies. Only a few results will be presented here. 

6.1 Selection of h — > bb and Measurement of M{ul) — M(xT) 

This analysis is very similar to the one for LHC Point 5, but since the signal cross sections 
are smaller, harder cuts are needed: 

• $ T > 100 GeV; 

• > 4 jets with p T > 50 GeV, p T ,i > 250 GeV p T ,2 > 150 GeV; 

• Transverse sphericity St > 0.2; 

• M cS > 1000 GeV; 

• $ T > 0.2M eff . 



Events with exactly two tagged b jets were then selected. Figures [34] and [35] show the bb 
mass distributions for Points 1 and 2 respectively. The light Higgs would be discovered in 
this mode at either point with 10 fb~ or less, although h — > 77 would still be needed to 
provide the most precise mass determination. 

Events passing these cuts and having a bb mass within 2a of the peak were then selected. 
These events were also required to have two and only two additional jets with p T > 100 GeV. 
This cut is quite inefficient, but attempts to improve its efficiency produced more background. 
Each of the two jets was combined with the bb pair. The smaller of the two masses was 



selected for each event and plotted in Figures.|3q and |3y]. The endpoints for the decay 
sequence qi — > X2Q ~^ Xi^bq are 739 GeV and 751 GeV respectively. Measurement of this 
edge will be limited by statistics to ~ 50 GeV for 10fb _1 . 

6.2 Selection of %i ~^ Xj^ + ^~ 

Sleptons are quite heavy for LHC Points 1 and 2, so there is no edge in the dilepton mass 
spectrum as for Point 5. There is, however, aZ-» £ + £~ signal that can be used to distinguish 
these otherwise rather similar points. 

The basic selection cuts given at the beginning of the previous subsection were applied. 
In addition, the events were required to have two opposite-sign, same-flavor leptons with 
Pt > 10 GeV, \r}\ < 2.5, and E T < 10 GeV in a cone R = 0.2. The mass of the two highest 



Pt such leptons is shown in Figure |38] and Figure |39] for LHC Points 1 and 2 respectively. 
The Z peak is questionable for Point 1 but rather clear for Point 2, which has larger tan/3 
and hence more mixing of gauginos and Higgsinos. Even for Point 2, however, the statistical 
error on the number of Z's is ~ 15%. Clearly this is a measurement that needs higher 
luminosity. 
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There are also observable signals in the like-sign and in the opposite-sign, opposite-flavor 
dilepton channels. These would certainly be useful in a global fit, but they do not seem to 
provide the sort of precise measurements being considered in this paper. 
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Figure 34: M(bb) for pairs of b jets for the LHC Point 1 signal (open histogram) and for 
the sum of all backgrounds (shaded histogram) after cuts described in the text. The smooth 
curve is a Gaussian plus quadratic fit to the signal. The light Higgs mass is 111.2 GeV. 
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Figure 35: Same as Figure E£|for LHC Point 2. The light Higgs mass is 125.1 GeV. 
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Figure 36: The smaller of the two bbj masses for the signal and background events with M(bb) 
within 2a of the peak in Figure |34| and exactly two additional jets j with px > 100 GeV. 
The endpoint of this distribution should be approximately M^ ax = 739 GeV. 
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Figure 37: Same as Figure |3<| for LHC Point 2. The endpoint should be M^ ax = 751 GeV 
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LHC Point 1 




Figure 38: M« for the Point 1 signal (open histogram) and the sum of all backgrounds 
(shaded histogram). 
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Figure 39: M« for the Point 2 signal (open histogram) and the sum of all backgrounds 
(shaded histogram). 
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7 Determining SUSY parameters 



Once a number of quantities have been measured, we can attempt to determine the particular 
SUSY model and the values of the parameters. The strategy will be to attempt to perform 
a global fit to the model parameters using all of the available data, much as the Standard 
Model is tested using the W and Z masses and the many quantities precisely measured by 
LEP/SLC. Such a fit is beyond the scope of our work, and we adopt a simpler procedure. We 
assume that from measurements of global parameters such as those discussed in section |^ 
we know the approximate scale of the superpartner masses and have some idea that we 
might be in a SUGRA model. The object is then to determine the parameters of that model 
and check its consistency. We must therefore determine the parameters mo, m^, A , tan/3 
and sgn/i. As we will see A is difficult to determine. Its value is given at the unification 
scale, and the value that is relevant for the phenomenology is the one evolved down to the 
electroweak scale. Many choices of Aq evolve to the same (fixed point) value, so there is 
reduced sensitivity to its value at the unification scale. In addition it always appears scaled 
by the Yukawa coupling for the relevant quark or lepton. Hence its effect on lepton masses 
and on the quarks of the first two generations is very small. Its effect is significant for top 
squarks and at large tan (5 for bottom squarks. | 

Our strategy for determining the parameters is as follows. We choose a point randomly in 
parameter space and compute the spectrum. We assign a probability to this point determined 
from how well it agrees with our "measured quantities" using our estimates of the errors on 
those quantities. The process is repeated for many points and the probabilities used to 
determine the central values of the parameters, their errors and their correlations. The 
probability distribution functions are not always Gaussian; the + ler (— lcr) errors quoted 
below are such that 15.87% of this distribution is below (above) the quoted value. Thus, 
68.27% of the probability falls within our definition of ±ler. 

7.1 LHC Point 3 

At LHC Point 3 (mo = 200 GeV, m 1( / 2 = 100 GeV), LEP will discover h and measure its 
mass. The relevant error is that from theoretical calculations of the mass in the supergravity 
model which is likely to dominate the error from the LEP measurement. We will assume 
an error of ±3 GeV. Using the results presented in Section [| we assume the following 
measurements and errors: 

• M^o - M^o = 52.36 ± 0.05 GeV, 

• Mg-M- b = 20.3 ± 2.0 GeV, 

• M h = 68.3 ± 3 GeV. 

As described above the mass difference M~ g — M~ h is insensitive to the mass assumed for Xv 
Using the strategy outlined above we get the following constraint on the parameters: 

*Th.e Higgs mass that we use is that determined from 1-loop calculations as implemented in IS A JET [Q. 
Two loop corrections to the Higgs mass, implemented, for example, in SPTHYIA |l3| lower the mass by 
about 5%. 
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• mo = 200±l 3 GeV, 

• mi/2 = 99.9 ± 0.7 GeV, 

• tan/? = 1.99 ±0.05, 

• the sign of \i is determined to be —1, 

• A Q is constrained to be greater than —400 GeV. 

There are no clear correlations between the parameters. The additional constraint that the 
average value of the light squark mass is within ±20 GeV provides no additional restriction 
on the parameters. The relevant phenomenological parameter is the value of A evolved down 
to the the electroweak scale. The relevant values are those for the third generation, A^ and 
A t . Information on these can only come directly from data on bottom and top squarks. In 
this case A t = — 176 ± 22 GeV. Ab is not well constrained; it is allowed to range from 50 to 
— 500 GeV. Over this range the mass of the two stop eigenstates varies only slightly from 
270 and 320 GeV at one end to 260 and 330 GeV at the other. Constraining A^ is very 
challenging as it appears scaled by the small factor of the bottom quark Yukawa coupling. 

The degree of precision may be surprising. Over most of the SUGRA parameter space 
and Xi are gauginos (i.e., they have no Higgsino components), M^o — M^o then determines 
rrii/2- Mg and M^o are then predicted and a consistency check of the model made by the 
measurement of M g . Information on mo, is then obtained from M~ h and the value of is 
sufficient to constrain tan/3. 

The other measurements available at this point are now used to provide powerful con- 
sistency checks of the model. The measurement of the bottom squark and other squark 
masses and the event rates for isolated leptons without jets discussed in Section |3| are ex- 



amples. Another example is the branching ratio for \\ ~^ Xi^ + ^ • Figure ^0] shows this 
branching fraction for a selection of SUGRA models that have parameters in the slightly 
larger range m = 200 ± 15 GeV, m 1/2 = 100. ± 1.5 GeV, tan/5 = 2.0 ± 0.1. At Point 3 
BR(x2 — * Xi e+e ~) — 16.5%. This branching ratio can be constrained using the method 



described in section 3.4;. 



7.2 LHC Point 5 

At LHC Point 5_(m = 100 GeV, m 1/2 = 300 GeV, tan/3 = 2.1), h will be discovered at LHC 
in its decay to bb from its production in the decays of supersymmetric particles, and its mass 
will be measured precisely from its decay to 77. Using the results presented in Section ||, we 
assume the following set of measurements: 



M^l - Mf/M^l - M*o/Mf = 108.6 ± 1 GeV, 
The decay g — > tt is allowed, 



The end point of the spectrum in Figure g3| is 506 ± 40 GeV, 
M h = 104.15 ±3 GeV. 



44 



These results correspond to two possible solutions: 



m = 100.5^ 2 GeV, 
mi /a = 298^ 6 GeV, 
tan/3 = 1.8±g;i, 

A* = ±i; 



and 



• m = 91 ± 3 GeV, 

• m V 2 = 288 ± 18 GeV, 

• tan (3 = 3.1 ±0.2, 

• /i = —1. 

Both of these solutions provide good fits to the "data" . While there is no constraint on A, 
the values of A t and A b are constrained: A b = -740 ± 180 (A b = -750 ± 220) GeV and 
A t = —495 ± 30 (A t = —536 ± 56 ) GeV for [i positive (negative). The main differences 
between the mass spectra for these two solutions are in the masses of x%i X% X%i an d xt 
and the masses of the heavier Higgs bosons; they are significantly larger in the positive /x 
case. The claimed sensitivity to the lepton decay spectrum (see Figure ^) of 5m ~ 5 GeV 
implies that this parameter's range can be narrowed somewhat. We have not investigated 
the sensitivity of this spectrum to the sign of fi. At this point the errors used (except for 
the one on Mh) are limited by statistics, so additional luminosity will cause the errors to 
drop. Using 0.6 GeV for the error on the dilepton endpoint and 23 GeV for the error on 
the endpoint of Figure ^3| reduces the errors on m 1( / 2 and m to ±7 and ±9 (±2.5 and ±10) 
respectively for positive (negative) \i\ the error on tan (3 is not reduced. 

As well as the differences in the masses of the heavier Higgs and gauginos noted above, 
the negative sign solution has larger stop masses and hence a smaller branching ratio for 



g — ► tit. We have investigated the sensitivity of the top quark signal discussed in section |5/4 
to the sign of \x. The parameters m = 90.4 GeV, mi/ 2 = 290 GeV, tan/3 = 3.1, \i = —1, and 
A = 66 GeV have the opposite sign of fi from Point 5 but give an excellent fit to the "data" 
used for fitting at this point.| A sample of 100K events was generated for this parameter set, 
and the analysis leading to Figure p3| was repeated. The result is shown in Figure HTt note 



that the bins in this figure are twice as wide as those in Figure |33|. It can be clearly seen that 
the amount of reconstructed top in SUSY events is reduced relative to that at Point 5 and 
that this fact can be used to eliminate this alternative solution. In addition, the solutions 
with fj, = — 1 have smaller branching ratios of q — > qx\h so the observed number of higgs 
events should be able to severely constrain this case. 



§This modified point has a combined probability of 98% of fitting the data. 
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7.3 LHC Point 4 



At LHC Point 4 (m = 800 GeV, m 1/2 = 200 GeV, tan/? = 10), determination of the 
parameters cannot be done by the simple method already described. Here only two masses 
can be measured in a straightforward manner; the light Higgs mass from its decay to 77 and 
the Xi ~ X2 mass difference from the endpoint in the dilepton mass distribution. We use 
these two measurements and the determination of the SUSY scale from the M e g analysis of 
Section ||: 

• M^o - M=o = 69 ± 1, 

• M h = 117 A ± 3 GeV, 

• mm(M~ g , M Qr ) = 580 ± 60 GeV from the M eff analysis. 
These constraints restrict the parameter space to two regions: 

. m = 784±^ GeV, 

• mi/a = 200 ± 8 GeV, 

• tan (3 = 9 ± 2, 

• A is not constrained, 

• A* = +1, 

• A b < 200 GeV and A t < -150 GeV; 



or 

• tan/? = 14 ±4, 

• m = 950 ± 210, 

• mi/2 = 185 ± 10 GeV, 

• A is not constrained, 

• jj, — — 1 

• A b = -160 ± 150 GeV and A t = -400 ± lOOGeV. 

The uncertainty on M^o — M^o is limited by statistics at low luminosity. A reduction in its 
error would reduce the error on m^. 

In order to constrain parameters further, models whose parameters are consistent with 
these values would need to be generated and their predictions for the distributions shown 



in Figures 17, 20, and Ell calculated. Those which are inconsistent with the "observed" 



distribution can then be rejected. This exercise is beyond the scope of this paper. 
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7.4 LHC Points 1 and 2 



At LHC Point 1 and 2 (m = 400 GeV, m 1/2 = 400 GeV, tan/? = 2, 10), event rates are low 
and precision measurements difficult at low luminosity. We use the following constraints. 

M h = 111.4 ± 3 GeV for Point 1, or 

M h = 125.4 ± 3 GeV for Point 2, 

min(Mg, M Qr ) = 920 ± 90 GeV from the M eff , 



and 



The end point of the spectrum in Figure £3| is 745 ± 50 GeV. 
In the case of Point 1, there are two solutions: 
m unconstrained, 
m 1/2 = 400lf8 GeV, 
tan/3 = 2.0±g;f, 
fi = +1, 

A is not constrained, 

A b = -1100 ± 200 GeV and A t = -650 ± 55GeV; 
one with negative /i that has instead 
m 1/2 = 392±fg GeV, 
tan/3 = 3.3±g;|. 

In the case of Point 2, there is again a solution for either sign of fi: 
mo unconstrained, 
m 1/2 = 405^? GeV, 
tan/3 = 10.6 ±0.3, 
/x = ±l, 

A is not constrained, 

A b = -1100 ± 200 GeV and A t = -800 ± 55GeV. 

The lack of a constraint on mo at these points is alarming but can be explained, mi/2 is 
large and the renormalization group scaling from the GUT scale forces the squark masses to 
be comparable to the gluino mass almost independent of the input value of ttiq. In cases of this 
type one needs to measure slepton masses which are less affected since their renormalization 
group scaling is controlled by a wea k rather than a s . Since the slepton masses are of order 
500 GeV, this is a difficult task. Even a lower bound on the masses would constrain itiq. 
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If we reduce the error on the end point of the spectrum in Figure |23] to 20 GeV which 
might be achievable with high luminosity running at LHC, the uncertainty on ra\M reduces 
to ±30 GeV. As in the case of Point 5, the decay g — > tt is allowed. An analysis similar to 
that discussed there should be able to establish that this channel is open. 

In the case of these points, the negative \l solution can be eliminated. The branching 
ratio for \\ Xi% is shown in Figure [42] as a function of tan f3 for both signs of \x for the 
solutions in the allowed range. It can be seen from this Figure that the branching ratio is 
substantially larger for ji = —1. The decay chain — > qx® —> q^Xi leads to a small Z 
peak as shown in Figure |38| which corresponds to BR{\2 ~^ Xi^) = 0.6%. The smallest 
branching ratio for the \x = — 1 solution shown in Figure f|2| is 5.7%, approximately a factor 
of ten larger. It is clear that these two cases can be distinguished in 10fb _1 despite the fact 
thata larger data set may be required to measure the size of the peak shown in Figure 38. 
In the case of Point 2, The difference in branching ratios is even greater, again enabling 
elimination of the \i — — 1 solution. 
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J, Point 3 B(x 2 -> Xi e + e ) 




0.1 0.2 0.3 

B(X 2 -> Xi e + e") 

Figure 40: Branching ratio x° — * Xi^ + ^~ f° r the SUGRA models m = 200 ±15 GeV, 
m 1/2 = 100 ± 1.5 GeV, tan/? = 2.0 ± 0.1. 
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Figure 41: The sideband subtracted invariant mass distribution of Wb candidates. The 
analysis is as described in section |5]3] The solid histogram is for Point 5 and is the same as 
that shown in Figure ^ except that the binning has been changed. The hatched histogram 
corresponds to the point m = 90.4 GeV, mi/s = 290 GeV, tan/3 = 3.1, sgn/i = —1 and 
A = 66 GeV. 
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Figure 42: The branching ratio X° -> X?^ for models with M h = 111.5± 3 GeV, M eff = 
920 ±90 GeV and both signs of /i = +1 (circles) and /i = — 1 (crosses) as a function of tan/5. 



50 



8 Conclusions 



In this paper we have outlined a strategy that can be used to systematically explore super- 
symmetry assuming that it is discovered at the LHC We have given an example of a global 
variable (M e g) that can be used to determine the mass scale of SUSY if nature has chosen 
the SUGRA model. Such global variables will be used to give the first indication of a signal 
independent of the type of SUSY model. The production of heavy mass states virtually 
guarantees that events with very energetic jets will exist. Other SUSY models such as those 
with R— parity violation []TJ] may have no missing E? if the LSP decays within the detector. 
If the LSP decays to leptons, then all SUSY events will have leptons and they can be used in 
a global variable. If the LSP decays to hadrons, the jet multiplicity will rise and a variable 
similar to M e g should be effective. 

More detailed exploration will depend on the particular SUSY model. In the SUGRA 
models, there is a general feature, namely that the second lightest neutralino almost always 
decays to h + Xi if the channel is open and to £ + £~ + Xi with a substantial branching 
fraction if it is not. In the former case, this will be the dominant source of h and it will 
be discovered in this process via its decay to bb if it has not been seen at LEP. In the 
latter case the measurement of the position of the end point in the £ + £~ mass distribution 
provides a very precise measurement of the mass difference between two of the sparticles. 
After first observing one of these signals, one will move up the decay chain to determine 
other quantities. 

We have then illustrated, using specific examples, some techniques that can be used to 
determine masses and branching ratios of sparticles. Some of these quantities were then 
used to determine the fundamental parameters of the SUGRA model some of which can be 
determined with great precision. The ultimate goal of such studies would be to use very 
many measurements to make an overconstrained fit to the model, rather in the same way 



that current data are used to test the Standard Model fT5 

The results in this paper are only an indication of the exciting physics that lies ahead for 
the members of the CMS and ATLAS collaborations if nature proves to be supersymmetric 
on the weak scale. 
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